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Abstract 

We investigate small— x resummation effects in QCD coefficient functions for 
Zog and Wg fusion processes, and we compare them with the known ones of 75 
type. We find a strong process dependence, that we argue to be due to the possible 
presence of collinear singularities for either small or large k of the exchanged gluon. 
For top quark production, we find that the gg tt and Z^g — > tt channels have 
larger resummation enhancements than the Wg ^ tt one. 
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1 Introduction 



In the past few years, various applications have been developed [1-5] of the k-factorization 
method to combine high-energy behaviour in QCD with the renormalization group 
{RG). 

On one hand, for small— a; hard processes of strongly interacting partons, resummation 
formulae have been derived for various coefficient functions and for the next-to-leading 
quark entries of the singlet anomalous dimension matrix. Resummation effects turn out 
to enhance scaling violations in the HERA energy range ^, and could be responsible 
for the small— X rise of structure functions seen experimentally P]. 

On the other hand, k-factorization is also relevant for the spontaneously broken elec- 
troweak theory |Q, where it provides a generalization of the effective W approximation 
for high-energy weak boson fusion processes, like, e.g., Zqq —>■ ti and Wg —>■ tb. In 
such cases, it takes into account correctly important high-energy off-shell effects which 
are neglected in the (naive) collinear approximation [ITU . 



The purpose of the present note is to calculate in detail the coefficient functions of 
Zqq and Wg fusion type by concentrating on gluon off-shell effects, and to compare the 
corresponding resummation formulae with the one of 7(7 type previously found. Note that 
while for a virtual photon source both transverse and longitudinal coefficient functions are 
available, here we limit ourselves to longitudinal gauge bosons, whose exchange dominates 



the fusion process, and is related - by the off-shell equivalence theorem [^, 11| - to the 
corresponding Goldstone boson exchange. 

It will turn out eventually that resummation effects are strongly process dependent, 
and do not enhance much the Wg tb channel compared to, say, the Zog, gg ti ones 
for top production. Part of the paper is devoted to explore the reasons for such a fact, 
and also to provide a quick way to estimate enhancement factors in the various cases. 

Roughly speaking, off-shell effects yield a cross-section increase because the phase 
space for the exchanged gluon transverse momentum k opens up at high-energies, and 
is furthermore weighted by the anomalous dimension exponent 77V — (Sas/vr) (A^ — 
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N being the moment index. Therefore, for large enough energy s, the integration region 
where |k| exceeds the hard scale Q [Q"^ ^ ^ s) becomes increasingly important, if 
not suppressed by the squared matrix element. 

Let us refer to the phase space region Q < |k| as "the disordered k" region, to 
distinguish it from the normally ordered region Qo ^ |k| < Q, typical of the incoming 
parton jet. We find that the size of the cross-section enhancement at high energies is 
related to the possible existence of a collinear singularity in the disordered k region, 
where the exchanged quark aligns to the hard probe, rather than to the parton in the 
opposite direction. We also give the splitting functions of GLAP type which describe 
quantitatively the hard off-shell cross section, thus generalizing what already noticed in 
the photopro duct ion case [|I|. 

For instance, while the DIS hard cross section = a {'jig Qq) - due to longitu- 
dinal photons - is not particularly enhanced, the one for Z^g —>■ ti is instead normally 
enhanced, because of the collinear behaviour of the Goldstone boson QQ process. The 
corresponding splitting function yields a reasonable estimate of the enhancement factor. 

On the other hand, a {WlQ — ^ ib^ is not particularly enhanced, despite the existence of 
similar collinear properties of the ly— Goldstone boson. But in this case, since mf, ^ nit, 
the small k-region is also enhanced by a log {mt/mf,) factor due to the customary collinear 
singularity in the parton jet. Thus the ratio of large k to small k effects becomes of order 
unity. 

The contents of the paper are as follows. In Sec. 2 we set up the calculation of the weak 
boson-gluon coefficient functions on the basis of our previous treatment of k-factorization 
in electroweak processes. The actual analytical computation is performed in Sec. 3, and 
in Sec. 4 we compare the present results with the ones for the photon source. Here we 
also discuss the role of the disordered k region, and we provide the relevant splitting 
functions in the various cases. Our results are discussed in Sec. 5, and some analytical 
details are left to Appendices A and B, where we also discuss the relevant longitudinal 
splitting functions. 
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2 k- Factorization for V-g fusion processes 



We consider here contributions to heavy quark production in which the fusion subprocess 
is of the mixed type, e.g. Z^g ti, Wg tb, the hard scale being essentially the top 
quark mass (Fig.l). Since we have already discussed off-shell effects in the weak- 
boson channel for on-shell gluons, we concentrate here on the strong interaction effects, 
and the correspondig resummation formulae of QCD corrections. 

The fusion process that we consider (Fig. 1) is, to start with, described by double- 
k± factorization, in which the hard subprocess cross section a (v (q) g (k) — > (51Q2) is 
factorized from the structure functions in the hadron (lepton) at hand. This means that 
in the high energy regime 

s^mj^,q^,k^, —,— fixed, (2.1) 
s s 

the Q1Q2 production cross section takes the form 

^QiQ2 ^i^) = J d'^q.dyd'^kdzTv (y, q) J^g (z, k) avg f-i, ^ j . (2.2) 



Here J-'y (J^g) are the unintegrated V{g) structure functions of the initial particles, and a 
is defined by a projection with eikonal vertices ( appropriate in the high energy limit ) of 
the off-shell Vg Q1Q2 absorptive part A^^, in a physical gauge A ■ n = 0, as follows 



In the following, we assume J-'v to be given by the lowest order V-emission probability 
off the external fermion, i.e., 

jr ) = 9tt9A \^ (1 _ y) (2.4) 

On the other hand, the unintegrated gluon density, J^g {z, k) will also contain the QCD 
radiative corrections that we want to take into account, and is defined in terms of the 
gluon-hadron absorptive part in the gauge p ■ A = in the usual way [Q 

„ U2 

k) = / -—.—P'-fG,. (k, p) (2.5) 
(^27rj zs 



so that 

/"^ d^kJ^g{z,\^)=g^{z,Q^) (2.6) 

•J 

is the corresponding gluon density in the initial hadron A. 

The k-factorized expression (2.2) shows off-shell effects due to the dependence of a on 
the transverse momenta q (k) of the weak boson (gluon). The q-dependence was investi- 
gated in Ref. and amounts to calculable corrections to the effective-W approximation 
in the high energy limit. The effect is not too large, and was estimated to be of the order 
of 10% for the present value of the top mass. 

On the other hand, the k-dependence is enhanced for s ^ m^, or z <^ 1, by large 
anomalous dimension effects in J^g. They come from QCD perturbative contributions of 
the form z~^gn (k"^, Qo^) (a^log (l/z))'^, which are resummed by the BFKL equation |T^ , 
and are translated in the coefficient function by the k-integration in Eq.(2.2). 

In order to concentrate on strong interaction effects, we shall take in Eq.(2.2) the 
"small q" limit for the weak boson, so as to translate it into a single-k factorization 
formula. This is done by setting, for q = 0{My), 

~ yp^ + q^ ^ yp^ (2.7) 



and by using in Eq.(2.3) the Ward identity |Q, |TT 



2/ W^Sp<. ^ ^'^"^l^^p. = Mv'A<;^, (2.8) 

which relates the longitudinal V amplitudes to the ones for the corresponding Goldstone 
bosons, and yields 

q Wv9\ci^=0 = ^^TfT^^pa = — , TT^ (2.9) 



2ys |k2 

We then notice that the q^ integration in Eq.(2.2) can be explicitly done around |g| ~ Mi 
by using Eq.(2.8) and the expression (2.4) for jFy, to obtain 



(J, 



We have thus exhibited the 0(1/My-) part of the cross section (2.2), due to the 
longitudinal V polarization, that will be considered in the following, and is also dominant 
in the region \q^\ = O (My), provided My -C m^. The additional O (l/m^) terms can in 
principle be fully evalueted from (2.2) on the basis of the perturbative analysis of the full 
douh\e-k± cross section (2.3), but are known not to provide large off-shell effects 0] and 
will no longer be considered in the following. 

The single-k factorized formula (2.10) can be further elaborated to provide a resum- 
mation formula for the corresponding coefficient function. 

Firstly we define the Mellin transform in the energy variable 

(2-11) 

ij J 

to rewrite Eq.(2.10) in the form 



X r°° ds 


s 


' Jo s 


(mi + m2)^. 



where the 1/ (A^ — 1) factor comes from the y integration, and we have defined M = 
nil + 1712- 

Secondly, we express the k-dependence of jF/v in terms of the BFKL anomalous di- 



mension 



1 /k^^ 7]v(«3) 



as follows 

(k) = (as) (^^j 4 M , (2.14) 

where /i = O (M^) is the factorizaton scale. 

Finally, inserting Eq.(2.14) into Eq.(2.12) allows performing the k-integration in terms 
of the calculable k^-moments 



k2 \M2/ \M2'm2 
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and provides the final result 

J^^h^"^' {^N {as {M'))) ■ gi (M^) . (2.16) 



-Q1Q2 f2\ _ 1 / 9v + 9a 



Here the factor in front of gf^ provides the QCD coefficient function for the given pro- 
cess. Since •jn in Eq(2.13) is a known function of the effective coupling as/ (A^ — 1), the 
expression (2.15), evalueted at {o^s {M'^)), provides the resummation formula we are 
looking for, once the lowest order expression for the off-shell cross section o"*^^ is given. 

Note that, in deriving Eq.(2.16), we have kept, for symplicity, frozen at its value 
as (M^), and we have also used the expression (2.14) even for < Q^, where Qo is a 
scale defining the boundary of the perturbation approach (a^ (Qq) < 1). It can be proven 
|l|, |13[ however, that using a RG improved expression (2.14) and/or a full solution of the 



BFKL equation [12] for jFjv(/c), including higher twists and possibly running coupling 
jl^ , does not change the final result (2.16), except for subleading terms of relative order 
as {M"^), which are not considered here. 

Thus our procedure will be to first evaluate the perturbative expression for the off-shell 
a^^ and then the corresponding /i— function, providing the resummed coefficient. 



3 Evaluating the coefficient function 

The hard sub-process under study — Vlq — > Q1Q2 — is a variant of the customary 
DIS process and has the simplifying feature of having only abelian-type diagrams (Fig. 
2). While the corresponding matrix element is thus rather straightforword, the analytical 
final state integrations are non trivial due to the unequal mass kinematics and will be 
sketched in the following. 

In order to calculate a'^^ in Eq. (2.9) it is convenient to introduce the off-shell trans- 
verse gluon polarizations 

= ^ (k'^ + ^P') ' 'out = (0, 0), (3.1) 
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satisfying e ■ p = e ■ \i. = 0, and the longitudinal one 



satisfying Sa ■ k = 0, and 

£a-£/3 = riJap, VaS^e^fi = g^y (3.3) 

with rjif^ = rjouT = —Vl = — 1- By using the Sudakov parametrization 

/ b;2 \ 1,2 

it is easy to realize that, in the high energy limit, 

e1-e'^^ = ^^lf. (3.5) 
Therefore, the eikonal projection in Eq. (2.9) becomes (with u = 2q ■ p) 

after noticing that Ai ijq = Auy ^ = by the usual invariant decomposition [0, ^ 

The absorptive parts occurring in the r.h.s of Eq. (3.6) are in turn obtained from the 
amplitudes M{G{q) +g{k) Qi{Pi, cri) + Q2{P2, (^2)) by the customary integrations over 
phase space and spins 

All + AjN,iN = E / dHW^ + iMj^^. (3.7) 



o"l,cr2 



The squared amplitudes in Fig. 2, summed over spins, turn out to be 



2 

V /I /\/f |2 , I \2\ 



Y^i\ML\' + \M^ 



IN\ 



\[m2 — t)[mi — u) J 
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{ml — m\Y 




+ 



(3.8) 



where Cv{Ca) denote the vector (axial) V Q1Q2 couphng constants, and we have 
defined u = 2q ■ p and the Mandelstam variables of the hard subprocess 



In the particular case of Zqq tt, the heavy quark masses are equal (mi = 7712 = m) 
and all terms in Eq. (3.8) drop out but the first, which yields 



Therefore, there is no explicit k— dependence of in this simple case. 

On the other hand, in the case of Wg tb, which has larger cross section for top 
production |]l5l, all the terms in Eq. (3.8) contribute (with mi = rrit and m2 = m^, say) 
and the k— dependence is more involved. 

In order to deal with the general mass configuration, it is convenient to write the two 
body phase space in terms of rescaled invariants, as follows (Appendix A) 



s = ik + qy, i={q-P,f = {k-P2f, 



u={q-P,f = {k-Prf. 



(3.9) 




(3.10) 



^$(1,2) 



(2^ 



5+((g + A)2 



ml)5+{{k-^f 





(3.11) 



where we have performed a trivial azimuthal integration, and we have defined 



m'? — i 



1 - r 




(3.12) 



T = 



ZV 



zv 
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It is also helpful to introduce the energy type variable 

a = r(l — t)zi' — ml{l — r) — mlr > ^^k^) 
e = 4r(l - r) 

in terms of which the phase-space boundary takes a mass-independent form: 



(3.13) 



d$(l,2) 



,6(06(1 -e)e(4a-ek^ 



(3.14) 



167^^/^^ 

It turns out that, by eliminating zu, t and u in favour of cr and ^, Eq. (3.8) takes, 
after some algebra, a particularly simple form 



2zv 4 



{C\{mi + m2f + Cy(mi - ms)^) 



(cr mf)(cr m\) 



+ 



.^9 ^9^'C. 9 9^9(-2(J + k2(|£- 1)) 



(3.15) 



'4^--^ ---^ ((7 + m?)2((7 + mi)2_ 

This expression has the remarkable property that the cr— dependence essentially factorizes 
from the ^—dependence, so that the various integrations needed are decoupled, eventually. 

For instance, the = 1 moment of acg, relevant for the high energy limit of Eq. 
(2.12) becomes 



-9l 



1 rf^ 



i67ryr^ya5 

4 



da 



C\{mi + m2f + C^(mi - ms)^ 
{a + m\){a -\-ml) 



((7 mi)2((7 m^)^ 



(3.16) 



The corresponding k^— moment, i.e., the /i— function in Eq. (2.15) is easily performed 
by doing the (linear) k^— integration before the cr and ^—integrations, with the result: 

m2\ 2ci;<j /mi -|- m2\ -27 1 



X 



(CKmi + ms)' + (mi - 7712)') 



] -5(1/2, 1-7)5(1-7, 1 + 7) X 

27 2 
mi — mn 



3-27 



(mi^ + 



27 
ma' 



2 r/;f 





+ 


- ml 






K) 


- Ill.y 


- mi ^ 



(3.17) 
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The equal mass case, relevant for Z^g —>■ QQ, is again particularly simple because only 
the first term in square brackets contributes. We obtain the mass-independent expression 



4ii(7, 1) 



asTR 2r(i-7)2r(i + 7)v/^ 

-ttG 



An 



r(3/2 - 7) 



(3.18) 



Of particular interest for estimating the size of scaling violations is the behaviour of 



Eq. (3.16) around 7 = (collinear limit) and 7 = 1/2 (saturating value [|T6| of the 
BFKL anomalous dimension). The 7 = limit is strongly dependent on the mass ratio 
r = 7712/1711 and has the form 



3 \ 1 _ fj 



(C^-Cl)logi 



1 — r 



(3.19) 



On the other hand, the large 7 behaviour is roughly determined by the presence of a 
7 = 1 double pole, with behaviour 



uG9 1^ 



TT 



1 — r\ 2 



1 + r 



Y^2 



(1-7)^ 



(3.20) 



to which, in Eq. (3.17), only the first term in square brackets contributes. The saturating 
value at 7 = 1/2 is instead given by 



TT 



4ii(l/2,r) = ya. 



1 — r\2 



1 + r 



(3.21) 



to be compared with the rough estimate obtained by setting 7 = 1/2 in Eq. (3.20). 

The overall behaviour of hP^[^^ is plotted in Fig. 3 for various mass ratios. In the 
equal mass case it increases from the 7 = limit all the way up to the double pole at 
7 = 1. The enhancement ratio at 7 = 1/2 (the asymptotic value at the BFKL Pomeron 
singularity) is 



r=l 



TT 



2.47, 



(3.22) 
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a value typical of other single— k processes, like heavy flavour photoproduction or DIS 
scaling violations [0, |]. 

Thus we conclude that the gZ —>■ QQ fusion process has rather large resummation ef- 
fects, while being disfavoured at Born level with respect to the gW and gg fusion channels 
for top production. 

On the other hand the unequal mass gW tb case is peculiar because the 7 = 
coefficient (Born cross section) is enhanced, according to Eq. (3.19), by a large factor 
~ logmt/mb due, as we shall see in Sec. 4, to a collinear singularity in the rrib = limit. 
Thus the /;,— function starts decreasing away from 7 = 0, has a minimum, and then is 
driven up to larger values by the 7 = 1 double pole (Fig. 3b). 

As a consequence, the enhancement ratio 



is not large, and finite energy resummation effects are not likely to be important. 

The above discussion shows that the relative importance of the Zog vs Wg channel 
for top production is energy dependent, the latter being dominant at low energies, but 
less enhanced by resummation effects. In any case, the most important process for top 
production is expected to be gluon-gluon fusion which yields a greater cross section at 
low energies and is expected to be strongly enhanced by QCD resummation effects at 
energies of LHC type (Cf. Sec. 4C). 

4 Process dependence of resummation effects 

We have already noticed that resummation effects due to large small— x anomalous di- 
mensions are substantially different in the Zog ti vs Wg — > th cases. Comparing with 
the known coefficient functions of (heavy) quark production with transverse |]l| and lon- 
gitudinal [0 photon sources, we find the enhancement ratios h{l/2)/h{{)) listed in Table 
1. 




(r 



30 



1 



(3.23) 
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We can thus roughly distinguish "normally enhanced" processes [ZlQ, JtO) from 
not enhanced ones {Wlq, 7^(7). We would like to understand this fact on the basis 
of the observation |l| that coUinear singularities in the internal momenta may enhance 
the k^— moments of the hard probe cross section. 

Since we compare large 7 with small 7 moment indices, we should distinguish at least 
two cases. 

A) Small k enhancement. 

This occours in some massless limit, e.g. rrig — > in DIS-type processes, or mh -C rrit 
in the cases examined here. The collinear singularity is due to a (nearly) massless quark 
exchange of virtuality i in the region k^ <^ \i\ <^ Q^, where Q (or M) is the hard scale. 
It yields the logmt/nij, factor in Eqs. (3.19) and (3.23), and a logQ^/k^ factor in the 
DIS case, where it provides a I/7 pole of the coefficient function /i2(7) (Cf. Ref. [3]). This 
enhancement, due to a "normal" collinear behaviour in the target jet, is important in 
order to assess the magnitude of the nearly on-shell (or 7 = 0) region. 

B) Disordered— k enhancements 

This is the typical off-shell effect that we are exploring in the high-energy regime, 
in which the k— phase space opens up away from the normal collinear region mentioned 
before. 

We shall call "disordered k" region the one in which Q^(M^) -C |t| <^ k^ < zz/ and 
therefore the exchanged quark aligns to the electroweak boson rather than to the incoming 
hadron. In such region the off shell cross section ang for a hard source {V or 7) coupled 
to quarks, is again dominated by a collinear singularity as follows 

Here, however, r = k.'^/zu is the Bjorken variable of the exchanged quark, as probed 
by the hard gluon, and Pu^qq is the corresponding splitting function. 
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The expression (4.1) provides the correct double-pole residue at 7 = 1 of the various 
coefficient functions in Table 1, where also the corresponding definition of the coupling 
gjj is quoted. In fact, the N — 1 moment of (4.1) is 



(|)=^.|^.|log^, (k^»g^), (4.2) 



Ah= ['dfPH^,^{T), (4.3) 

Jo 

and thus its contribution to the /i— function (2.14) becomes 

V=i(7)7^i ^ (i_^)2 ' (4-4) 

in agreement with the properly normalized expressions quoted in Table 1. 

Here we notice that the source dependence comes from both the coupling constant g\i 
and Au^ the N — 1 moment of the splitting function. It is easy to reahze (Appendix B) 
that for transverse photon coefficient functions we have (Table 2) 

P,,^,,(f)=f2 + (l-f)^ (4.5) 
while for longitudinal weak bosons (Goldstone bosons) we have 

PG^,,{f) = 1, (4.6) 

as also quoted in Table 2. 

It is clear that the double pole approximation provides most of the enhancement factor 
only in those cases in which the 7 = region is not enhanced. This explains the "normal" 
enhancement of 7^5' and Zlq processes and the small anhancement of WlQ processes. 

Particular attention is needed for the 7^(7 process, occourring in the DIS longitudinal 
structure function, which is substantially different from the corresponding case in broken 
gauge theories. In fact, there is no Goldstone boson in this case, and the region (A) is 
not coUinear singular due to the helicity flip zero at = for Breit frame scattering 
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of massless spin 1/2 quarks. However, the disordered region (B) happens to have a 
hnear coUinear divergence, being essentially of Coulomb scattering type (Appendix B). 
Its contribution to (Tl (defined by saturating with el in Eq. (B.5)) is of type 



aa 



= — 4rl-r, r = —, 4.7 

TT ZV ZV 

and provides a single pole behaviour of the /i— function 



KM ^ (4.8) 

TT i — 7 

Once again, the coUinear analysis sketched above explains the weak enhancement in 
this case, even if in a slightly more involved way. 

C) s— Channel enhancement 

Finally we should recall [jl[, for completeness, that non abelian gluon-gluon fusion 
processes show a third coUinear region, the one in the s— channel, which is important to 
estimate the size of resummation effects, e.g. in hadroproduction of heavy flavours. This 
is a phase space region of the hard subprocess gg QQ in which <^ s <^ zv and 
<^ P = (q + k)^ ^ ~ q^. So that both k^ and q^ are large and of the same order. 
Therefore, the massless quark limit is relevant, the cross section Cgg behaves as 1/s, and 
its flrst moment is dominated by the intermediate gluon coUinear singularity 

^--^[jp)=^^^'^Jp (4-9) 
for one produced flavour. Furthermore, this time 

Ah = J drPg^ggir) (4.10) 
is related to the gluon — quark pair splitting function in the flnal state. 
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It turns out that, at extreme energies, in which both gluon transverse momenta q 
and k carry a large anomalous dimension 'Jn{(^s), the behaviour (4.9) causes a triple pole 
in the coefficient function 



which occours at precisely the asymptotic value 7 = 1/2 because of the double— depen- 
dence. 

Of course, this implies that the non-abelian s — channel region is asymptotically 
dominant, so that the gg — » QQ channel is more strongly enhanced than the 'yg and Vg 
ones. 

5 Discussion 

There are various outcomes of the analysis presented in this paper. The ffist one is that 
the relative importance of the Wg vs the Zog fusion processes for top production changes 
with energy, when scaling violations drive the coefficient away from the 7 = region, 
collinear singular for Wg tb only. However, the gg —* ti process is even more enhanced 
by scaling violations, and is thus conffimed as the most important at energies of LHC 
type. 

The second point is that the enhancement factors are roughly described by the dis- 
ordered— k collinear singularities, which are process dependent, since they involve splitting 
functions of the hard source (transverse or longitudinal 7's, Z's and W^s in the present 
case). This leads to a variety of enhancement factors, and implies that it is not really 
possible to eliminate all large resummation effects by a redefinition of the gluon density. 

This observation is relevant for the analysis of small x scaling violations being per- 
formed at HERA. While it is fruitful to look for a factorization scheme -or definition of 
the gluon density- which incorporates some universal enhancement factors, one has to 




(4.11) 
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live with the fact that resummation effects are sizeable for some processes, and may in 
fact be needed to explain the small x rise of structure functions. 

Thus, a comparative study of various processes (e.g. F2,Fl, heavy flavour production) 
is needed at both theoretical and experimental level. 

There is a final point to be noticed. The coUinear analysis performed in this paper, and 
compared to exact results, provides in fact a method which can be generalized to gluon- 
gluon kernels and other processes for which exact results are not yet available. In fact, 
it provides a quick way to estimating resummation effects and their process dependence, 
even for non abelian hard subprocesses. This analysis is left to future investigations. 
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Appendix A: Unequal mass kinematics 

For the process G{q) + g{k) Q1Q2 we use a Sudakov parametrization of momenta, 
which, in the notation of Fig. 1, reads 



= yp' 




(A.2) 



A'^ = {-q + Pi)'' ^{k- P2Y = ZTp^ - yfp" + A' 



(A.3) 



where y{z) is the momentum fraction of the probe (gluon) with respect to the incoming 
lepton momentum p (hadron momentum p) . 
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Then the mass-shell conditions Pf = ml yield 

.(1-^)^'"' + ''^ + '^'', (l-.),= !!i±(^^ (A.4) 

zv zv 

where v — 2q- p — ys. Therefore, the expression of the invariants in terms of longitudinal 
and transverse variables is easily obtained. The i/— phase space is described by 

+ (A + q)^ ^ mj + (k - A)^ ^ + (A + q)^ ^ mj + (k - A)^ 

T 1 — T 1 — f f 

The first of these expressions is important for the normal coUinear kinematics, the second 
one for the disordered— k region. The remaining invariants are 



i^{q- P,Y = m? - zi/T + q2 + 2q • A = _ ^i^ + (^ + ^q)' _ ^2- (^.6) 

1 — r 

u = {q- P2f = ml- zv{l - r) + q2 + 2q • (k - A) (A.7) 

j^^^_(k + q)2>!!i + J!^. (A.8) 

T 1 — T 

Eqs. (A. 6) and (A.7) explain the definition (3.12) of r (for q^ = 0). By replacing 
(A.4) in the two-body phase space, and by performing the f integration first, we obtain 

^' ' r(l - r) 2(27r)2 V ^ i _ ^ ^(i _ ^) ) ^ ' 

and by then performing the d^A integration we arrive at the expression (3.11) of the text. 
Finally Eq. (A. 6) provides the phase space boundary 

|£| = |A|2>ff^ + q2) (A.IO) 
^ 1 — r ' 

relevant for the disordered— k kinematics, in which, by (A. 5) f ~ \s? /zv. 
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Appendix B: Longitudinal and scalar splitting func- 
tions 



Here we want to derive the splitting functions Ph qq relevant for the coUinear behaviour 
in the disordered— k region, because they may involve longitudinal and scalar polarizations 
(possibly in a broken gauge theory), which are not usually treated in the literature. 

With the notation of Fig. 1 and referring to the splitting process V{q) Qi{q + A) 
+(52 (—A), we consider the probe momentum q^^ to be off-shell and we use a Sudakov 
parametrization with massless quarks 

q'^^P^-^jf, {q^^-Q\2p-p^v), (B.l) 

-Pf = (g + A)^ = (1 - f)p^ + jT^^p' + A^ {{q + A)2 = 0), (B.2) 

A 2 A 2 

p'M ^ = fp'^ - ^ ^-A^ i^^^-z ). (B.3) 

(1 — t)v ^ 1 — 

where we have reabsorbed the momentum fraction y in p'^, for simplicity. 

Here the relevant region of "disordered— k" is m^, -CA^ -C k^, so that we will 
eventually also set Q^/k^ —>■ 0. 

The vector boson polarizations in the Landau gauge parallel the ones in Eqs. (3.1)- 
(3.3) for the gluon and are given by 

e^*, = (0, 0), = (0, Bout: 0) (B.4) 

with £ouT • A = 0, and by 

£^: = ^(?'^ + 2^p'^), (el^l). (B.5) 

Correspondingly, the quark (antiquark) spinors with helicity A(— A') take the usual form 



, AV(-A) , 

v^{-A) = ^fE,\ J \ , (B.6) 

77^ (-A) 
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(B.7) 



where the relevant bispinors will be approximated, in the coUinear region by the expres- 



and similar ones for with 9 = |A|/|g(l — r). Note that off-shell effects are only relevant 
at order 6"^. 

Since we only consider, for the various initial polarizations, the leading coUinear singu- 
larity, the splitting functions can be defined according to the probabilistic interpretation, 
by summing over quark helicities in the usual way: 



sions 





A 





2A2 



A,A' 



i:iM(£^^AA')r 
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Transverse polarizations 

In this well known case there is no much difference between broken and unbroken 
theories. The massless quark limit is smooth, the lowest order matrix elements, given in 
Table 3 

Mxx' = u\q + A)^{Cv + CAl^y'i-A) (B.IO) 

are helicity conserving and show the well known zero in the forword direction, due to a 
clash with angular momentum conservation. Correspondingly, the definition (B.9) has a 
finite = limit, yielding the customary splitting function 

PT^,,{f)^f' + {l-ff, g'^^C'y + Cl (B.ll) 

as quoted in Table 2. This one corresponds to the normal logarithmic coUinear singularity. 
Longitudinal and scalar polarizations 

Here we should single out the broken theory case, in which the q'^ term in the longi- 
tudinal polarization (B.5) dominates the vector exchange contribution, for small , as 
explained in Sec. 2. 

By picking up the contribution, restoring quark masses and using the Ward Identity 

M^y ~ ^u^ [Cv{mi - 1712) + CA{mi + 7712)75] v^' (B.12) 

we end up with effective (pseudo) scalar couplings, which are singular in the Q ^ limit 
and helicity violating. 

The 1/Q factor in (B.12) is reabsorbed in the definition of the Goldstone boson cross 
section aog in Eqs. (2.8) and (2.9) , while for the (pseudo) scalar matrix element we can 
use the massless quark kinematics as before, thus obtaining the (pseudo) scalar values 
of Table 3. Note that there is again a zero in the forword direction, due to a clash with 
angular momentum conservation, which in this case would imply conserved helicity (!). 
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As a consequence, in the Goldstone case the coUinear singularity is again logarithmic, 
with constant splitting function 



i^G^qq[T) - 1, Qg - Cs + Cp - La > (B.idj 

as quoted in Table 1. Note that in the definition of a^^ in Eq. (2.9) and of h in Eq. (2.16) 
an overall factor of (mi + 1712)"^ /My has been taken out. 

The longitudinal case in the unbroken theory is peculiar also. In fact, this time the 
qi^ term in (B.5) yields vanishing contribution because of current conservation, and the 
p'^ term is of order Q, but is helicity conserving. As a consequence, there is no clash with 
angular momentum conservation and no zero in the forword direction. By using Table 3 
and Eq. (B.9) we then obtain 

Pw = 4f2(l-r)^, \gL\'^C'y + Cl (B.14) 

As a consequence, the fact that the longitudinal process has to vanish on-shell (Q^ = 
0), is compensated by a linear coUinear divergence as stated in Eq. (4.6) of the text. 
Integration over it, taking into account the phase space boundary (A. 10), provides an 
effective splitting function 4f(l — f) but no logarithm and thus, upon f— integration the 
same residue 2/3, but a single pole, in /i(7).This result is connected with the well-known 
fact that longitudinal polarizations yield subleading logs in the coUinear analysis. 

To summarize, the longitudinal polarizations yield: (a) the (scalar) helicity flipping 
contribution (B.12), due to current non conservation, which in broken gauge theories yields 
the dominant contribution, the probability density being of order (mass)^/(5^| A|^; and 
(b) the helicity conserving contribution (B.15), the only one contributing in the unbroken 
case, whose probability density is of order 1/|A|^. The corresponding splitting functions 
are given in Table 2. 
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Table 1 



process 


t^9h 


P{z) 


"' \ 1 J 

coupling 


It 

HO) 


h(l/2) 

MO) 


Zg^QQ 


Ml 


1 


(yJI'n 1 

-ft 

TT (1-7)2 


2.00 


2.47 


Wg^tb 


^[(l+r)2Cl+(l-r)2c2] 
w 


1 


TT (1 — 7)2 


1.33 


1.49 


irg QQ 


2 


^2 + (1 - 


a;i«2 1 

TT 3(1-7)2 


1.71 


2.48 


It QQ 


Tre^ 


z^ + {\- zf 


as'-/kW/2 1 

TT 3(1-7)2 


4.00 


4.00 


ll QQ 


Tre^ 


4z(l - ^) 


1 

TT 3 1 — 7 


1.00 


1.45 



Table 2 



i 


P^(z) 


\coupling constant\^ 


S 


1 


ci + cl 


Ll 


4^(1 - z) 


cl + cl 


IN 


(1 - 2z)2 


ci + c\ 


OUT 


1 





^Helicity conserving part only. The helicity flipping contribution is equal to the scalar 
contribution with C5 = (mi — m2)/Q and Cp — (mi + m-2)lQ. 
2 Effective splitting function after phase space integration (see App. B). 
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Table 3 



A 


A' 


i 


^ ab 


+ 




S 


k 

-{Cs-Cp) 




+ 


S 


{Cs + Cp) 1 


+ 


+ 


L 


-2{Cv + CA)Q^Jz{l-z) 






L 


2(Cy - CA)g\/^(l - z) 


+ 


+ 


IN 


-{Cv + Ca)^=^={1-2z) 
Jz{l - z) 






IN 


T 

{Cv-Ca) , ^ {l-2z) 

^Jz{l - Z) 


+ 


+ 


OUT 


-i{Cv + Ca) I ^ ^ 
^z{l - z) 






OUT 


i{Cv Ca) I 
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Figure Captions: 



Fig. 1: Kinematics of tlie Vg contribution to heavy quark production. 



Fig. 2: The lowest order amphtudes for the hard sub-process gV^QiQ2. 



Fig. 3: The function h{'y) for (a) mi = m2 and (b) 1112/ rrii ~0.03, relevant for single 
top production. 
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Fig. 2 
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